Scattering of dark-state polaritons in optical lattices and quantum phase gate for 

photons 



M. Masalas^'^ and M. Fleischhauer^ 

^ Technische Universitat Kaiserslautern, 67663 Kaiserslautern, Germany 
Vilnius University Research Institute of Theoretical Physics and Astronomy, 2600 Vilnius, Lithuania 

(Dated: February 9, 2008) 

We discuss the quasi 1-D scattering of two counter-propagating, dark-state polaritons (DSP), each 
containing a single excitation. DSPs are formed from photons in media with electromagnetically 
induced transparency and are associated with ultra-slow group velocities. State-dependent elastic 
collisions of atoms at the same lattice site lead to a nonlinear interaction. It is shown that the 
scattering process in a deep optical lattice filled by cold atoms generates a large and homogeneous 
conditional phase shift between two individual polaritons. The latter has potential applications for 
a photonic phase gate. The quasi 1-D scattering problem is solved analytically and the infiuence of 
degrading processes such as dephasing due to collisions with ground-state atoms is discussed. 
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A major challenge for quantum information processing 
using individual photons as qubits is the implementation 
of logic operations. Such operations require efficient non- 
linear interactions for pairs of photons which cannot be 
achieved in conventional optical materials. As the ef- 
fect of the photon-photon coupHng depends on the non- 
linear susceptibilities as well as the interaction time, it 
has been suggested to use ultra-slow light in resonant 
systems with electromagnetically induced transparency 
(EIT), where the interaction time is long and the non- 
linear susceptibilities become large due to resonance en- 
hancement 0, 0, Q. At the boundary of an stationary 
EIT medium Hght pulses become spatially compressed 
in the propagation direction by the ratio of group ve- 
locity "Ugr to the vacuum speed of light As a con- 
sequence the number of photons in the pulse decreases 
by the same factor. Excitations are temporarily trans- 
ferred to the medium by the formation of quasi-particles, 
so-called dark-state polaritons (DSP) which are superpo- 
sitions of electromagnetic and atomic degrees of freedom 
Q. As another consequence of the pulse compression the 
interaction time in a head-on collision stays constant irre- 
spective of the value of Ugr. Thus in order to achieve long 
interaction times co-propagating pulses were considered 
0. In this case the interaction is however not homoge- 
neous and it is difficult to avoid spectral broadening of 
the wave packet. 

We here suggest a completely different mechanism 
for an efficient nonlinear interaction between ultra-slow 
light pulses. The slow-down corresponds to a shift of 
the polariton composition from pure photons to matter- 
waves Furthermore the pulse compression leads to 
an increasing density of the matter component. Thus 
collisional interactions between atomic excitations can 
yield an effective nonlinear coupling between two wave- 
packets. To further enhance the strength of this in- 
teraction we consider a lattice potential in the tight- 
binding limit. Starting from a fully quantized effective 
one-dimensional model of Hght propagation in a lattice 
we derive analytic solutions for the quantum scattering 



of two single-photon wave-packets in the s-wave scatter- 
ing limit. We show that the pulses attain a homogeneous 
conditional phase shift which may be large enough for 
the implementation of a quantum phase gate. 




FIG. 1: top: Atomic five-level system with quantized probe 
fields E± of opposite circular polarization and propagation 
direction. Q, denotes the Rabi-frequency of the classical, un- 
depleted control field, bottom: The atoms are assumed to be 
confined in a 3-D lattice potential with lattice constant a. In 
the tight-binding regime they occupy only the lowest Wannier 
state of effective width Zq. 



Let us consider a cold gas of bosonic, five-level atoms 
as shown in Fig. in a deep three-dimensional lattice po- 
tential under tight-binding conditions. The atoms form 
an M-type system, with the ground state \g) and the 
excited states |e±) coupled by two orthogonal polariza- 
tions of a quantized probe field E± propagating in the -\-z 
or —z direction respectively. The excited states are fur- 
thermore coupled to meta-stable state \q±) by a classical 
probe field of Rabi-frequency Vl. All atoms are initially 
in the ground state \g). The atoms are described by five 
Bose fields -i/'i, where i £ {g,e.±,q±\ denotes the inter- 
nal state. The Hamiltonian of the system is given by 
H = i/at + -ffiat + ^^at-f + ^^coU, whcrc the atomic part 
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reads 
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+ nuj,]ij^. (1) 



-ffiat is the lattice potential, which is assumed to be the 
same for all internal states. The interaction with e.m. 
field reads in rotating-wave-approximation 



ydV ^ Vl.?if^(r)V'9. +H.C. (2) 



Finally we take into account collisions between atoms 
in the internal states \g) and \q±) as well as coUisions 
between atoms in \q^) and \q+) in s-wave approximation 
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CoUisions of atoms in the excited state are irrelevant as 
these states will attain vanishing small population. Here 
£'£^'' (r) are the positive frequency parts of the probe field 
operators corresponding to the two orthogonal polarized 
modes, ft (r) is the Rabi frequency of the control field, 
and the it's describe the collision strength, which can be 
expressed in terms of the corresponding s-wave scatter- 
ing length Ui = ATraih/m. The Wi's are the frequencies 
corresponding to the electronic energy levels. 

Using eqs.QJ-lEll, and ijSJ the Heisenberg equations for 
the field operators can easily be obtained: 



= (^1^^' + v + nLOg)i^g- ;u*i?i-Ve 
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n(ugki^l^Pg + u±kijl^'ipqAiij. 



(4c) 



We here have not included decay from the excited states 
as these states will attain only negligible population. 
The above operator equations are nonlinear and thus im- 
possible to solve exactly. Therefore approximations are 
needed. 



First it is assumed that both probe and control fields 
propagate along the z axis and can be written in the fol- 
lowing form: E^^'^ {r,t) £:±(r,i)e'(±'=^-"*),17(r,t) = 
no{t)e-"^''\ Here both £±ir,t) and ^oit) are slowly 
varying functions of r and t. One also assumes that the 
probe field coupling ^£±/Tl is much weaker than the con- 
trol field coupling Qq. In this Hmit we may consider the 
control field undepleted and classical and can thus set Qq 
constant. The assumption that all atoms were initially 
in the ground state g then also implies that only a small 
fraction of atoms is excited to the states q± . 

The strength of the lattice is considered to be large 
enough such that only the lowest energy level of each 
potential well is occupied and tunneling between the 
wells is negligible. In this limit the atomic field op- 
erators can be expanded in the basis of (real) Wan- 
nier functions Q. Since only the lowest lattice state 
is occupied, only Wannier states of the lowest Bloch 
band Wj{r) = VFo(r — r^) survive in this expan- 
sion. In a deep lattice Wannier-functions of neigh- 
boring sites have only negligible overlap and one has 
J d'^rWi{r)Wk{r) — Sik. This expansion along with 
a separation of fast oscillating terms yields ipg = 

and = T,j ^3 (i") g±j (i) e±*('=-'=<=)^-*("-"'=)* . Here 
the operators g^, e±j, and q±j are slowly varying in t 
and the summation runs over all lattice sites. 

Substituting the Wannier-expansion into eqs.l^ yields 
equations for the slowly varying atomic variables. Since 
the Wannier functions are well localized, the slowly vary- 
ing functions remain almost constant within one lattice 
site. 

5j = (i-j) e-j + ij^X (rj) e+j- 

(ws5].9j + Wg+^lj-g+j + ^h-^-fl-i) -^9j (5a) 



e±j = -iAe±e±j + i^£± (tj) gj + i^oq±j (5b) 



q±j = -iAq±q±j + in'^e±j (5c) 

-i("±±9±jg±j + + ug±g]gjj -^q±j 



(4b) with the detunings Ae± = uJe± 
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LUq± — u! + uJc + fi{k — kc) /2m. We here have set LOg = 0. 
The factor / = a/lo describes the confinement strength 
in the lattice and is defined as / d^r Wi(r)^ = P /a^. 

In order to solve the above equations of motion for the 
matter-field operators a weak-probe approximation will 
be applied. In zeroth order of the probe field the states 
|e±) and \q±) remain unpopulated and one finds for the 
ground state operators 
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(6) 
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If the lattice has a regular filling with a well defined 
number Nj = N oi atoms per site, we can make the 
replacement gjgj — > in eq.Q. A regular filling can be 
achieved e.g. by employing a Mott-insulator transition 
in a lattice |a,0|- this case the self-phase modulation 
described by ij^J can just be absorbed in the definition 
of the energy of the ground state atoms. Furthermore 
we can replace the ground-state operator in the equa- 
tions for the operators of the other states by a c-number 
gj ^/N . With this we obtain in first order of the probe 
field 



■A -u-^^^^- 
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are strongly localized around the center of the poten- 
tial wells, the microscopic polarization changes rapidly in 
space. On the another hand in Maxwell's equations for 
the electric field only the macroscopic polarization enters. 
The macroscopic polarization can be obtained by aver- 
aging over a volume small compared to the wavelength. 
If the lattice constant a is sufficiently smaller than the 
relevant wavelength of the probe field, the lattice struc- 
ture disappears in the polarization. In the opposite limit 
effects Hke Bragg scattering can occur, but we are not in- 
terested in these phenomena here and thus assume a < X. 

Using the macroscopic polarization we find the follow- 
ing equation of motion for the slowly varying amplitude 
of the field operator 



q±] = -i(Ag± -I- Ug±f^n)q±j + inle±j 

-jw±±^9±j9±i9±j - iu+^-^ql-jq^jq±j (7b) 

where n = N/ is the average density of atoms . Next 
we assume resonance conditions, i.e. Ae± — and Ag± -|- 
Ug±np = and that the probe field varies sufficiently 
slowly. Under these conditions we can apply an adiabatic 
approximation. 

In zeroth order of the adiabatic approximation the time 
derivative in the equations of motion for the atomic vari- 
ables is neglected which yields 



(8a) 



with £±j = £± {vj). Thus 

(0) i+\ t j-3 W?n M c c 

-nu^-f^^£U^,£±, (9) 

Proceeding in the same manner with the first order of 
the adiabatic approximation we find that the operators 
corresponding to the excited states |e±) contain a term 
proportional to the time derivative of the probe fields 



it) = eg (t) 
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Here higher order terms in 8±j containing a time deriva- 
tive were neglected as they correspond to corrections 
higher order in both the probe field and the adiabaticity 
parameter. 

The adiabatic solutions for the matter fields can 
now be used to calculate the slowly varying am- 
plitude of the probe-field polarizations P±(r, i) = 
^l* X), I^j(i")p3±i(0e±i(0- Since the Wannier functions 
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with group velocity Vgr = 2c?i|riopeo/lMp'^'^ (assuming 
fgr <C c). The appearance of the group velocity in the 
denominator of (lllll suggest on first glance an diverging 
nonlinear interaction when the group velocity approaches 
zero. One should take into account however that due 
the pulse compression the total photon number A^ph '-^ 
J d^r is only a fraction Wgr/c of the input value. Ijll|l 
becomes much more transparent if it is translated into 
an equation of motion of the dark-state polaritons Q 

^'±(z,t) = cos d£±{z,t) - sm9VAipq±{z,t) (12) 



where £ — £ y^huj/2eoA, A being the cross section of the 
light beam, Vgr — ccos^ 9, and ipq± are the slowly- varying 
amplitudes of the matter fields.. In the adiabatic limit 
considered here the orthogonal quasi-particle, the bright- 
state polariton is not excited and thus £± — ^± cos 9. 
With this we find the following propagation equation in- 
side the medium 
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_^2a+^^3^^^^^^^ (13) 



where we have substituted the s-wave scattering length 
Gij and the recoil velocity Wroc = haj/mc. The polariton 
number densities n± = ^'l^(z)\E'±(z) undergo a sudden 
increase at the boundary of the medium since the elec- 
tric field is continuous there. Inside the medium they 
propagate form-stable with Vg^. 

Eg . (11311 leads to a self- and cross-phase modulation of 
the dark-state polaritons. If initially only one polariton 
is excited of each sort the self-phase modulation vanishes. 
To solve the quasi-lD scattering problem it is convenient 
to introduce the two-particle wave-function 



w{z,z',t) = (O|*+(z,i)*_(z',i)|0) 



(14) 
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where \4>) is the initial state vector of the system and |0) 
corresponds to the polariton vacuum. One can show that 
in the case considered here, namely only one polariton of 
each class is initially excited, all information is contained 
in w. In terms of center-of-mass and difference coordi- 
nates R = (z + z')/2 and — z — z' the equation of 
motion for w reads: 

This equation has a simple interpretation. The l.h.s. de- 
scribes the propagation of the two components in oppo- 
site directions. The r.h.s. describes an interaction for 
z — z\ i.e. when the two polaritons meet. The inter- 
action conserves the center-of-mass of the two polaritons 
and the solution of ifT^ll reads 

w{R,£_,t) = w(i?,^-2wgri,0)exp(-iA(/)e(O) (16) 

where &{£,) is the Heaviside step function. One recog- 
nizes that the shape of the two photon wave-function 
remains unchanged by the collisions and there is only a 
homogeneous collision phase. This is illustrated for the 
case A(j) = tt in FigEl 

The conditional phase shift between the polaritons 
which originate from two single-photon wave-packets is 
then transferred back to two photons at the exit of the 
medium. In this way a quantum phase gate between two 
individual photons could be implemented, if the condi- 
tional phase shift can reach the value of tt. The collision 
phase in itTfill is given by 



It is interesting to note that A(j> does not depend on the 
pulse parameters. It is not necessary that the two pulses 
have a certain length or the same shape. Both should 
however occupy the same transverse mode. One also rec- 
ognizes that the use of a lattice potential has two impor- 
tant effects. First of all phase diffusion of the individual 
dark-state polaritons caused by the scattering of atoms in 
states \q±) with ground-state atoms is ehminated by the 
regular filling. Secondly the local enhancement of the 
density leads to an enhancement factor f^. In a deep 
lattice / can be as large as 10. To give an estimate 
of achievable phase shifts let us assume — lOnm, 



Vg, = lOwrcc, A = X^, X = 800nm, and / = 10. This 
yields a phase shift on the order of unity, which is of the 
required order of magnitude. 

The main limitation of the present scheme is set by 
the finite dephasing time of the dark-state polaritons. 
For very small group velocities the interaction time T ~ 
L/2vgr of the pulses becomes rather long. Since the com- 
pressed pulse length in the medium L should be suf- 
ficiently larger than the wavelength, the minimum T 
is about lOA/wgr which for the above parameters with 
Vrec ~ 1 cm/sec gives T « 40/isec. In slow- light and 




FIG. 2: Snapshots of two-particle wave-function w{R,^,t) at 
equidistant times as function of center-of-mass coordinate R 
and difference coordinate ^ = z ~ z' . ^ < > 0) cor- 
responds to polaritons propagating toward (away from) each 
other. Overlapping components generate a phase flip A(j}, 
which is here taken to be n. 



light-"stopping" experiments 0, 0, dephasing times 
of msec have been observed and with the use of a deep 
lattice even longer times should be achievable. 

In summary we have shown that scattering of ultra- 
cold atoms in a deep 3-D lattice together with the trans- 
fer of excitations between photons and atomic excitations 
through dark-state polaritons can be used for a condi- 
tional homogeneous phase shift between individual pho- 
tons. An essential requirement to obtain sufficiently large 
phase shifts is the transverse focusing of the polaritons 
to a cross section comparable to and sufficiently long 
dephasing times. 
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